IP-ASTP-06-95 
CU-TP-679 
hep-th/ 9503200 



Quantum Nucleation of Vortex String Loops 

Hsien-chung Kacfi and Kimyeong Lee* 

^Institute of Physics, Academia Sinica 
Taipei, Taiwan 
and 

*Physics Department, Columbia University 
New York, New York 10027 



Abstract 

We investigate quantum nucleation of vortex string loops in the relativistic 
quantum field theory of a complex scalar field by using the Euclidean path in- 
tegral. Our initial metastable homogeneous field configurations carry a spacelike 
current. The path integral is dominated by the 0(3) symmetric bounce solution. 
The nucleation rate and the critical vortex loop size are obtained approximately. 
Gradually the initial current will be reduced to zero as the induced current inside 
vortex loops is opposite to the initial current. We also discuss a similar process in 
Maxwell-Higgs systems and possible physical implications. 
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We consider the theory of a complex scalar field with a global U(l) symme- 
try. When a field configuration is homogeneous, its conserved four-current can 
be timelike, spacelike or lightlike. Here, we consider quantum dynamics of clas- 
sically stable configurations carrying a spacelike current. We argue that they are 
quantum mechanically unstable, and decay through quantum nucleation of vortex 
string loops whose current reduces the existing current inside the ring, leading to 
the slow decay of the initial current. 

When there is a uniform external electric field, electron-positron pairs could be 
created via quantum tunneling and will fly apart because of electrostatic force [1]. 
In Yang-Mills Higgs systems with magnetic monopoles in broken phase, monopole- 
antimonopole pairs would be created in a uniform magnetic field [2]. If the gravity 
is taken into account, magnetically charged black hole pairs can be created when 
there is a uniform external magnetic field [3]. Similarly, in open string theories, 
where positive and negative charges are attached to the ends of the strings, open 
strings would be created by tunneling in a uniform electric field [4]. In the deSit- 
ter spacetime, there can also be a nucleation of closed vortex strings due to the 
expanding spacetime [5]. 

In the dual formulation of the theory of a complex scalar field in three-dimensional 
spacetime, vortices appear as charged particles and the uniform current does as a 
uniform external electric field. Thus, vortex-antivortex pairs will be nucleated and 
fly apart if there is a uniform spacelike current. The present work is to understand 
a similar process in four- dimensional spacetime. Instead of vortex-antivortex pairs, 
closed vortex string loops will be nucleated and expanded. In condensed matter 
physics, there is a long history of investigations on the vortex loop nucleation in 
superfluids in finite temperature due to the thermal and quantum effects, whose 
review can be found in Ref [6]. (Recent discussions about tunneling involving vor- 
tices can be found in Ref. [7].) The global charge current in superfluids is timelike 
and the thermal nucleation of vortices in superfluids is due to the relative velocity 
between superfluid and normal fluid. 
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In this paper, we consider the case where the initial current is spacelike. Clearly 
such a system cannot be realized naively in superfluids whose current is timelike. 
We, however, hope that this work enhances our understanding of the vortex string 
loop nucleation in general and other related issues like what configuration the 
system would settle down to after decay of the initial current. In a similar context, 
there has been a recent effort to understand the vortex string loop nucleation in a 
scalar field theory [8], whose analysis we feel is somewhat heuristic as the current 
the paper assumes is time-like. In addition, we also consider briefly the case when 
the scalar field is coupled to the gauge field. Additional physical implications will 
be discussed at the end. 

Our treatment of the vortex string nucleation follows the Euclidean path in- 
tegral method to calculate the imaginary correction to the initial energy density 
[9]. The vortex loop nucleation rate would then be twice the imaginary part of the 
false vacuum energy density. We will consider the leading order contribution and 
make some comments about the one-loop correction, which is much smaller than 
the uncertainty in our approximation to the leading order contribution. This work 
complements the work by one of us in which case the initial four-current is timelike 
[10], in which case the physics is completely different. The quantum nucleation of 
vortex ring is not possible energetically. If the initial configuration with time-like 
current or a uniform charge density is quantum mechanically unstable, it will decay 
by the bubble nucleation and expansion. 

This paper is organized as follows: We first discuss the characteristics of the 
initial configurations which is classically stable and spatially homogeneous. Then 
we discuss the Euclidean path integral including the vortex strings by using the 
dual formulation. We then solve the bounce equation approximately, which allows 
us to determine the radius of the string loop at the moment of the nucleation and 
the vortex nucleation rate per unit volume. We then investigate a similar case in 
a Maxwell-Higgs system and argue that the initial configuration with an external 
current could decay through the quantum nucleation of vortex string loops. Finally 
we conclude with some remarks. 
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A theory of a complex scalar field = /e*^/V2 is given by the Lagrangian, 
£=\{d»f? + \f(d^) 2 -U(f). (i) 

The conserved current for the global U(l) symmetry is J M = f 2 d^ip and its con- 
served charge Q = J d 3 xJ°. The ground state of the system is given by < / >= v, 
where the global symmetry is spontaneously broken. In this vacuum, there are 
massless Goldstone bosons and massive scalar bosons. In the broken phase there 
also exist topological global strings around each of which the phase ip changes by 
2tt. 

Here, we are interested in quantum decay of classically stable homogeneous 
configurations with a global current J M . Depending on whether the current is 
timelike or spacelike, we can choose a frame where there is a uniform charge density 
but no spatial current or a uniform spatial current but zero charge density. In 
this paper we consider the case when the current is spacelike. We then choose 
a reference frame where only nonzero component of is J z — j. For a given 
current j, the homogeneous f,ip configuration is given by / = u, ip — —jz/u 2 . 
This homogeneous configuration is static in time as we require the u to satisfy 

j 2 + u 3 U'(u) = 0. (2) 



To find out whether such a static homogeneous configuration is stable clas- 
sically, let us analyze the field equation under small fluctuations. Under small 
fluctuations around this configuration, f = u + 5f, ip = —jz/u 2 + Sip and 5/, 5i[) ~ 
gi(TOt-k-x)^ we g e ^. a dispersion relation from the field equation. The massive mode 
of mass m = a/j 2 /m 4 + U"{u) is always stable but the massless mode could be 
unstable. For small k, the dispersion relation for the massless mode is 

w 2 - k 2 + k 2 + ^»-3j 2 ^ ( ) 

w -k x + k y + u4u „ {u)+j 2 k z- W 
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This mode is stable only if 

u A U"{u) - 3j 2 > 0. 



(4) 



The two conditions, (2) and (4), can be obtained easily by considering an effective 
potential U e ff(f) = —j 2 /2/ 2 + £/(/). A stable homogeneous configuration appears 
as a stable point of this effective potential. It is intersting that the Goldstone mode 
with the z momentum behaves like a sound wave of speed less than unity. 

We now choose a specific potential 

U{f) = \{f-v 2 ) 2 . (5) 
For this specific potential (5), the condition (2) becomes 

j 2 + Xu 4 (u 2 - v 2 ) = 0. (6) 



When the magnitude of the initial current j is less than a critical current 
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. 2V\v 

3c = w (7) 



there are two solutions u± of Eq.(6), u- < a/2/3-u < u + < v. The stability 
condition (4) is satisfied only by the u+. The classically stable initial configuration 
we are considering here is then given by / = u + ,ip = jz/u\. When the initial 
current is very small, Eq.(6) implies that u + ps v(1 — j 2 /2\v 6 ) + 0(j 4 ). As the 
initial current gets larger, u + decreases to ^2/3v. There is no solution to Eq.(6) 
when \j\ > j c . Thus there is a maximum homogeneous current the system can 
endure. 

If the initial current is smaller than j c , a classically stable homogeneous config- 
uration is possible. However, we know it is always quantum mechanically unstable 
by quantum nucleation of vortex string loops whose interior current is opposite 
to the initial current. The standard method to calculate the decay of metastable 
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state by quantum tunneling is the bounce method in the Euclidean path integral. 
Once we know the bounce solution of the Euclidean field equation satisfying the 
correct boundary condition, the imaginary correction to the energy density of the 
metastable state is ImE/V = Ae~ B where the exponential suppressing factor B is 
the difference between the bounce action and the action for the metastable con- 
figuration and the prefactor is naively the ratio of the path integral prefactors. In 
the prefactor for the bounce solution, we take out the zero modes and take half of 
the absolute value of negative eigenvalue. The decay rate per unit volume is then 
T = 2lmE/V, where V is the volume factor. 

Since our initial configuration is expected to decay due to the nucleation of 
vortex rings and their expansion, it would be useful to describe the vortex dynamics 
more explicitly. This is possible in the dual formulation of the scalar theory [11]. 
For a given configuration of vortices, the ip(x) is multi- valued. The position of the 
a-th vortex would be given by qa(a a ) with the world sheet coordinates a a , a = 0, 1. 
For such a configuration, there is a conserved vortex current K^ v (x) = e pvp(J d p d ( jip, 
which can be expressed as 

= 2tt [ d 2 ae a Pd a q»dpq»5 4 (x p - q p {a)). (8) 

a J 

The orientation of the world sheet coordinates fixes the orientation of the vortex 
string, or the increasing direction of the phase variable ip. 

In the dual formulation, the generating functional in the original /, vp fields is 
rewritten in terms of the /, an antisymmetric tensor field C pv and the vortex string 
position qa(r,a). The on-shell relation between the original and dual variables is 
fd^ = \e» vpa H vpa where the field strength is H^ p = d^C vp + d v C Pil + d p C^. 
The dual Lagrangian is then given by 

A? = \W? + ^p H U ~ U(f) + \C^ V . (9) 

From the Lagrangian (9), the canonical momentum density for Cy is Hoij/f 2 . 
Since i?oi2// 2 for the homogeneous initial configuration is j /u 2 , the wavefunction 
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for the initial configuration is then proportional to exp \i J d^xjC^/u 2 }. In the 
dual formulation, the uniform current becomes the uniform "electric" field and 
exerts the "electrostatic" force on vortex strings, and so we expect vortex string 
loops to nucleate. 

The Euclidean path integral formalism is used here to calculate the imaginary 
shift of the energy density, which leads to the decay rate. We are calculating the 
Euclidean path integral, 

< F\e- HT \I >= J [f- 3 dfdC^dq^] V* F e- SE tt 7 (10) 

by the saddle point method. The Euclidean action can be obtained by the Wick 
rotation t = —it of the Minkowski action or by a direct dual transformation of the 
Euclidean path integral of the complex scalar field, and is given by 

Se = j d*x E |i(^/) 2 + 3p + U(f) - \C, V K^ . (11) 

Note that the interaction between vortices and C^ u is topological and is not affected 
by the Wick rotation. Since the initial and final configurations carry a nonzero 
current, there is a boundary term in the exponent of the path integral (10), 

£ = -i^p J d 3 x {C 12 (r F , x) - Ci 2 (r 7 , x)} . (12) 

Since the C^ v has the free boundary condition, the saddle point of the path integral 
will be the the stationary field configuration of the combined action, S = Se + £ 
which is also complex. From this combined action we get a bounce equation. The 
appropriate boundary conditions at r = ±oo is / — > u. 

Since the Euclidean initial configuration of the uniform current along the z axis 
has the symmetry under three-dimensional rotations in the r, x, y coordinates, we 
look for the bounce solution of this 0(3) symmetry. With r 2 = t 2 + x 2 + y 2 , the 
solutions would depend only on r, z. Since the boundary term is purely imaginary, 
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we see that the C^y field for the bounce solution should be pure imaginary to 
satisfy the boundary condition. We also choose the ansatz so that the C^y field 
has the spherical symmetry, 

C = \c^dx^dx v = iC{r, z) sin 9d9dip. (13) 

The field strength H = dC would be 



H = -H„ vp dx»dx u dx p = id r C sin 9drd9d<p + id z C sin 9dzd9dp . (14) 
6 



The boundary term (12) leads to the boundary condition C(r, z) — > jr 3 /3 for large 
r. For the field configuration to be smooth at the line r = 0, C(r, 2) should vanish 
when r = 0. 

From the translation symmetry along the z axis, we choose that the string lies 
on the z = plane and a two-dimensional sphere of a given radius R in the r, x, y 
space. The string current is invariant of the choice of the string internal coordinate. 
We choose the string coordinate so that q^(a a ) = (-Rsin a cos a 1 , R sin a sin a 1 , R cos a , 0). 
The delta function 5 4 (a^ - qP) becomes 5(r - R)5(6 - a )S(cp - a!)5(z)/R 2 sin 9. 
The string tensor current (9) becomes 

K^(x) = 27r(0><^ - § v <p^)S{r - R)S(z), (15) 

where 9^ = (9,0), cp^ = (<p,0). If it was the nucleation of a vortex with opposite 
vorticity along the z direction, there would have been a negative sign on the above 
vortex current. For the spherical symmetric solutions, the bounce field equation 
from the combined action is given as 



^d r (r 2 d r f) + d 2 J - -L {( drC f + (d z Cf) - U'(f) = 



(16) 



There cannot be a regular solution satisfying the boundary condition for arbitrary 
R. The reason is that we expect only one bounce solution, resulting in a unique 
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radius. Except for this value of R, we expect that the solution satisfying the 
boundary condition at the infinity will become singular along the line r = 0. The 
combined action for this configuration can be put as 

S E + Z = 4nJ drr 2 dz {^/) 2 + \{d z ff + ^RCf + (d z C) 2 ] + [/(/)} . 

(17) 

The bounce solution /fe(r, x), C&(r, x) of Eq.(16) can be analytically continued 
to lead to the solution of the Minkowski time. Let us first express the antisymmetric 
tensor field in the Cartesian coordinate, noting tan# = p/r with p = a/x 2 + y 2 , 
the two-form (13) becomes C& = iC^(r, z)(p 2 dr + pr/r 3 dp)dLp/r 3 . In the p > t 
region, the solution in Minkowski time is 

f(t,x) = fo(it,x) 

(2 
^dpchp - ^dtdif 

which leads to nonzero C p(fi , Ct^. By taking the exterior derivative of C, we can 
get the field strength. We can further analytically continue the above solution to 
the p < t region. From the Euclidean two-sphere covered by the vortex string, 
x 2 + y 2 + r 2 = R 2 , we get the radius for the vortex string on the x — y plane in 
Minkowski time as R(t) = y / ' R 2 + t 2 . The size of the string then increases with a 
constant acceleration. 

For a given initial current, it is not easy to find the bounce solution of Eq. 
(16) and get the tunneling rate. When the initial current is much smaller than the 
critical current, the size of the two-sphere would be much larger than the Compton 
wave-length of the Higgs particle and the / would be very close to u in this case 
away from the two-sphere. The bounce solution can then be well approximated by 
the solution of Eq.(16) with f — u outside the core region. Later we will show that 
this approximation is reasonably good. 

One conventional approach to this problem is to find a one-parameter A family 
of field configurations along the tunneling path and to calculate the effective energy 
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for these configurations, E(X) = K(\)\ 2 + V(X). We can then treat the problem 
by the one- dimensional WKB method. The tunneling path in the configuration 
space is then described by these one-parameter field configurations. In our case, 
we can have an effective string action coupled to the antisymmetric tensor fields. 
From this, we can get the energy functional for a vortex string ring as a function 
of its radius. However, the self-interaction term is complicated to treat. 

We follow another approach [2], which considers the one-parameter family of 
field configurations 0(A, x) in Euclidean time. Each field configuration is a solution 
of the Euclidean field equation, which might be mildly singular. Here, we imagine 
exploring the configuration space with these one-parameter stationary field config- 
urations. The Euclidean action S(X) of the configurations in this family is a finite 
function of this parameter and we find the stationary point among this family with 
respect to this parameter. The configuration at the stationary parameter would be 
the stationary configuration we looked for. The radius R of the two-sphere is the 
parameter. This approach is more profitable when there is a long-range force in 
the system as it treats the space and time components of the force simultaneously. 

When f = u, the bounce equation (16) becomes 

r 2 d r (J^d r C(r, z)^j + d 2 z C(r, z) = 2nu 2 R 2 5(r - R)5(z). (19) 

Since the equation is linear, we can require C to satisfy the boundary condition 
at infinity and to be zero on the line r = 0. We approach this by solving the 
homogeneous equation by the fourier transformation in the z coordinate, 

C(r, z) = ijr 3 + J dpe ipz (a p (pr - l)e pr + b p (pr + l) e - pr ) . (20) 

We impose the boundary conditions and match the homogeneous solutions at r = 
R. Our solution of Eq.(19) is then the sum of the background part 

Cback(r) = ijr 3 (21) 
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and the dynamical part 



C dyna = -^ j(it! 2 + r 2 + z 2 )m 



(R + rf + z 2 
(R-r) 2 + z 2 



ARr 



(22) 



There are several regions where we are interested in the detail behavior of Cd yna - 
Near the two-sphere, (R - r) 2 + z 2 « R 2 and 



C, 



u 2 R 2 



dyna 



In 



(R-r) 2 + z 2 

W 2 



+ 2 



(23) 



We can read off the delta function in Eq. (19) from this expression. Far away from 
the sphere, r 2 + z 2 » R 2 and Cdyna ~ — 4w 2 r 3 i? 3 /3(r 2 + z 2 ) 2 . For near the center 
of the sphere, r 2 + z 2 « R 2 and Cdyna ~ -4w 2 r 3 /3.R. 

Now we can ask how good our approximation f — u is. Far away from the 
two sphere, the C(r,z) approaches Cback and so we can see that / ~ u becomes 
the solution of the first of Eq.(16) due to Eq.(2). From Eq.(16) we can get the 
equation satisfied by the correction 8f — f — u due to the C^yna- Away from the 
two sphere, we can see Sf is approximated by 



2jd r C, 



dyna 



r 2 u 3 U"(u) 



(24) 



to the first order in C^ yna . We studied the asymptotic behavior of C& yna in various 
regions in the previous paragraph. We estimate 5f = —8j/(uU"(u)R) near the 
origin and 5f vanishes quickly for r 2 + z 2 » R 2 . Thus the correction is quite 
small if R is much larger than the mass parameter. We will later argue that the 
correction to the action due to 5f is also small compare with that from C(r,z), 
making the whole approximation consistent. 

As the / = u is a good approximation except near the two-sphere, our config- 
uration f = u,C = Cback + Cdyna 1S an approximately stationary configuration of 
the action away from the two-sphere. Because our approximate solution does not 
have any singularity at r = 0, we expect the singularity at r = of the solution of 
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Eq.(16) for arbitrary R » m will be mild. Our approximate solution is singu- 
lar near the two-sphere of distance less than m _1 . We now have a one-parameter 
family of the field configuration whose is approximately stationary except along 
the parameter R. 

Rather then calculating the action (17), which diverges at the infinity, we 
calculate the exponential suppression factor B, which is the difference between 
the bounce action (17) and the action for the initial configuration. We write the 
exponential factor B as a sum 



where H is the field strength for our approximate solution given by Eqs.(21) and 



The action B is the difference between the action for the approximate bounce 
solution and the background action. Since the C(x) is a sum of the background 
field Cback an d C'dyna; we get the B as a sum 



In the integration, we introduce the short distance cutoff e of order m near the 
two-sphere. This is because our approximate configuration becomes singular there. 
The first integral of Eq.(26) accounts for the interaction between the background 
'electric' field and the vortex string loop. The second integral is the self-interacting 
part of the vortex string loop. The integral can be performed by introducing two 
cutoffs (r — R) 2 + z 2 > e 2 and r 2 + z 2 < L 2 , by using the integration by parts and by 
using the asymptotic behaviors of the C- m t(x) fields in various regions. The short 
distance cutoff e near the two-sphere is of order mT 1 . Putting these together, we 




(25) 



(22). 




(26) 
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get 



B(R) = 47r 2 u 2 R 2 ln{cmR) - 2vr 2 j J R 3 , 



(27) 



where c is a constant of order one. 

Let us now consider what are the correction to the above equation. The solution 
of Eq.(16) is not singular near the two sphere and the correction to the factor B 
from the region inside the cutoff near the two sphere would be of order m 2 R 2 from 
the dimensional analysis. Since 5f ~ 1/-R, one can see the correction to B due 
to 5f would be again of order m 2 R 2 . Both corrections change the coefficient c by 
order 1, making Eq.(27) reasonable. 

Now we can find the stationary point of the exponential factor B with respect 
to R and get the bounce radius Rb » mT 1 as 

4-u 2 , / AJecmu 2 \ 
Rb « — In (28) 



3j V 3j 

and the exponential suppression factor at this radius becomes 

^^H ln HHr (29) 

By a simple scaling argument f/v , C/v 2 , ry/Xv , we can show that the classical 
action for the bounce is of the form Bb = jF(>j) with an unknown function F. 
Eq.(29) provides an approximate form for the function F. The radius (28) of the 
bounce solution, which is the critical radius, takes a form mR ~ (jc/j) m (j c /j) 
and the bounce action (29) is of order (j 2 /Aj 2 )(ln(j c /j)) 3 . Thus our semiclassical 
approximation is good if j << j c and A << 1. 

Let us now consider briefly the one-loop correction to our bounce solution. We 
can expand the fluctuations around the bounce solution in the spherical harmon- 
ics. The antisymmetric tensor field should be expanded in the tensor spherical 
harmonics. There will be four zero modes for the translation. There will be one 
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negative mode which corresponds to the change of the radius for two-sphere. There 
will also be gauge fixing terms for the antisymmetric tensor field. Since the orig- 
inal theory is renormalizable, we expect the dual formulation would somehow be 
renormalizable, even though the dual theory is a strong coupling theory. In prin- 
ciple, the ultraviolate divergence should cancel that from the background. From 
the standard technique [9,2], it is not hard to see that the normalization factor for 
each zero mode is of order 1/R. The one-loop correction would be then 

K = £ (30) 

with a dimensionless number a. The precise factor does not matter in our case 
because there is an uncertainty of order R in the B factor. By combining the 
results, the vortex string loop creation rate per unit volume becomes Ty = 2Ke~ Bb . 

In addition, we notice that a similar question can be addressed to the Maxwell- 
Higgs systems with the Lagrangian 

C = + \W) 2 + \f 2 ^ 9 + AJ 2 " U(f) - J^A,. (31) 

Without the external current, the system cannot be homogeneous in the presence 
of a uniform matter current due to the nonzero magnetic field. The initial config- 
uration is a homogeneous configuration f = u, A z = j /u 2 with an external current 
Jext = (0, 0, 0, j). The current carried by the Higgs field cancels the background 
electric current, resulting in the zero total current. For this homogeneous solution 
to be static the u satisfies Eq.(2). From the analysis of the small fluctuation around 
this configuration, one can show that this configuration is also classically stable if 
Eq.(4) is correct. As we are in the Higgs phase, there is no massless mode and in- 
stability can occur if the mass becomes negative. With the specific potential given 
in Eq.(5), such a homogeneous configuration can exist only if the external current 
is smaller than the critical current (7). The classically stable homogeneous config- 
urations can again decay quantum mechanically by the nucleation of local vortex 
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string loops. The analysis of such tunneling process would be somewhat similar 
to the global vortex nucleations in the previous discussions. We can use the dual 
formulation in the second paper of Ref. [11]. There is again an 0(3) symmetric 
bounce solution. The field configuration away from the two-sphere will fall expo- 
nentially to the initial configuration. Hence, there is no logarithmic term in the 
bounce action. The analog of Eq.(28) is B(R) = AnfxR 2 — 2iT 2 jR 3 , where // is the 
string tension in the small current limit. The bounce radius becomes = 4/i/3irj 
and the exponential suppressing factor becomes = 64/z 3 /(277rj 2 ). 

We have discussed the vortex loop nucleation in a uniform background current. 
We have analyzed the nucleation rate per unit volume and the size of the vortex 
loop at the moment of the nucleation. We argued that its size grows with a 
constant acceleration. We then extend our discussion to a Maxwell-Higgs theory 
with a uniform external current. We can approach our problem by an effective 
string action which contains the Nambu action and the interaction term between 
the background 'electric' field and strings. This would lead to a similar answer 
obtained previously if we take into account the string self-interaction. 

The two-fluid model for the superfluid works fine in describing the vortex loop 
nucleation in finite temperature [6]. Vortex loops nucleate thermally when the 
relative speed between two fluids is nonzero. The vortex nucleation in that case 
does not allow the naive path integral treatment as we do not know how to treat the 
normal fluid in the path integral. However, we can imagine a somewhat similar case 
where two superfluids exist at zero temperature, whose currents are both timelike 
but not identical. If we imagine a nontrivial coupling between them, vortices might 
nucleate. This process can be approached by the path integral. The characteristics 
of vortex loop nucleation will depend on the charge coupling. Since the current is 
timelike, there may be some surprise here. 

Finally, we can ask what is the final configuration of the string nucleations. (If 
we imagine a finite volume, string loops will be created and annihilate each other 
and disappear from the system. The initial current will be continuously reduced.) 
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Thus we expect the final configuration to be some sort of superfluid turbulence. 
In the theory with the local symmetry, the initial current of the Higgs field will 
decay away but the external current is frozen. We end up with a curious situation 
where there is the external current but no compensating current by the Higgs field. 
This configuration if static in time seems nonsensical because the static solution of 
the Maxwell equation has the the magnetic field growing linearly in space and so 
the energy density grows quadratically. The final configuration can be obtained by 
considering a similar system in two dimensional spacetime where the initial current 
decays by the bounce solution with vortex and antivortex lying on a Euclidean time 
axis. The initial Higgs current again decays away. The electric field grows linearly 
in time. However, one can show that the energy is conserved because there is no 
lower bound on the energy functional due to the external current. In short, there 
is no ground state in the local gauge theory when there is the external current. 
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note added: 

After this work appeared in a preprint form, we are informed by Dr. J.-M. 
Duan about his related works on the vortex nucleations in the charge density wave 
systems [12]. 
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